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ABSTRACT 

We have assembled a sample of high spatial resolution far-UV (Hubble Space Telescope Advanced 
Camera for Surveys Solar Blind Channel) and Ha (Maryland-Magellan Tunable Filter) imaging for 
15 cool core galaxy clusters. These data provide a detailed view of the thin, extended filaments in 
the cores of these clusters. Based on the ratio of the far-UV to Ha luminosity, the UV spectral 
energy distribution, and the far-UV and Ha morphology, we conclude that the warm, ionized gas in 
the cluster cores is photoionized by massive, young stars in all but a few (Abell 1991, Abell 2052, 
Abell 2580) systems. We show that the extended filaments, when considered separately, appear to be 
star- forming in the majority of cases, while the nuclei tend to have slightly lower far-UV luminosity 
for a given Ha luminosity, suggesting a harder ionization source or higher extinction. We observe a 
slight offset in the UV/Ha ratio from the expected value for continuous star formation which can be 
modeled by assuming intrinsic extinction by modest amounts of dust (E(B-V) ~ 0.2), or a top-heavy 
IMF in the extended filaments. The measured star formation rates vary from ^ 0.05 M© yr~^ in the 
nuclei of non-cooling systems, consistent with passive, red ellipticals, to ~ 5 M© yr~^ in systems with 
complex, extended, optical filaments. Comparing the estimates of the star formation rate based on 
UV, Ha and infrared luminosities to the spectroscopically-determined X-ray cooling rate suggests a 
star formation efficiency of 14lg^%. This value represents the time-averaged fraction, by mass, of gas 
cooling out of the intracluster medium which turns into stars, and agrees well with the global fraction 
of baryons in stars required by simulations to reproduce the stellar mass function for galaxies. This 
result provides a new constraint on the efficiency of star formation in accreting systems. 

Subject headings: galaxies: cooling flows - galaxies: clusters - galaxies: groups - galaxies: elliptical 
and lenticular, cD - galaxies: active - ISM: jets and outflows; stars: formation 



1. INTRODUCTION 

The high densities and low temperatures of the intr- 
acluster medium (hereafter ICM) in the cores of some 
galaxy clusters suggests that massive amounts (100-1000 
M0 yr~^) of cool gas should be deposited onto the cen- 
tral galaxy. The fact that this gas reservoir is not ob- 
served has been used as prime evidence for feedback- 
regulated cooling (see review by Fabian 1994). By in- 
voking feedback, either by active galactic nuclei (here- 
after AGN) (e.g., Guo et al. 2008; Rafferty et al. 2008; 
Conroy et al. 2008), mergers (e.g., Gomez et al. 2002; 
ZuHone 2010), conduction (e.g., Fabian et al. 2002; Voigt 
et al. 2004), or some other mechanism, theoretical models 
can greatly reduce the efficiency of ICM cooling, produc- 
ing a better match with what is observed in high resolu- 
tion X-ray grating spectra of cool cores (0-100 M© yr~^, 
Peterson et al. 2003). However, these modest cooling 
flows had remained unaccounted for at low temperatures 
until only recently. 

The presence of warm, ionized gas in the form of Ha 
emitting filaments has been observed in the cores of sev- 
eral coohng flow clusters to date (e.g., Hu et al. 1985, 
Heckman et al. 1989, Crawford et al. 1999, Jaffe 
et al. 2005, Hatch et al. 2007). More recently, it has 
been shown by McDonald et al. (2010, 2011; herafter 
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M+10 and M+11, respectively) that this emission is in- 
timately linked to the cooling ICM and may be the result 
of cooling instabilities. However, while it is possible that 
the warm gas may be a byproduct of ICM cooling, the 
source of ionization in this gas remains a mystery. A wide 
variety of ionization mechanisms are viable in the cores 
of clusters (see Crawford et al. 2005 for a review), the 
least exotic of which may be photoionization by massive, 
young stars. 

The identification of star-forming regions in cool core 
clusters has a rich history in the literature. Early on, it 
was noted by several groups that brightest cluster galax- 
ies (hereafter BCGs) in cool core clusters have higher 
star formation rates than non-cool core BCGs (John- 
stone et al. 1987; Romanishin et al. 1987; McNamara and 
O'ConneU 1989; Allen 1995; Cardiel et al. 1995). These 
studies all found evidence for significant amounts of star 
formation in cool cores, but the measured star forma- 
tion rates were orders of magnitudes smaller than the X- 
ray cooling rates (e.g. McNamara and O'Connell 1989). 
In recent history, two separate advances have brought 
these measurements closer together. First, as men- 
tioned earlier in this section, the X-ray spectroscopically- 
determined cooling rates are roughly an order of magni- 
tude lower than the classically-determined values based 
on the soft X-ray luminosity. Secondly, large surveys in 
the UV (e.g., Rafferty et al. 2006; Hicks et al. 2010), 
optical (e.g., Crawford et al. 1999; Edwards et al. 2007; 
BildfeU et al. 2008; McDonald et al. 2010), mid-IR (e.g., 
Hansen et al. 2000; Egami et al. 2006; Quillen et al. 2008; 
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TABLE 1 

Sample of 15 cooling flow clusters with MMTF YLa and HST FUV imaging 
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14h54m31.5s 


+18d38m32.4s 


0.0587 


0.025 


14.6 


39.0 


11980 


Abell 2029 


151il0m56.1s 


+05d44m41.8s 


0.0773 


0.040 


3.4 


527.8 


11980 


Abell 2052 


15hl6m44.5s 


+07d01ml8.2s 


0.0345 


0.037 


2.6 


5499.3 


11980 


Abell 2142 


15h58m20.0s 


+27dl4m00.4s 


0.0904 


0.044 


1.2 


< 2.5 


11980 


Abell 2151 


16h04m35.8s 
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0.0352 


0.043 


8.4 


2.4 


11980 


Abell 2580 ^ 


23h21m26.3s 
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46.4 
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Abell 2597 
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-12d07m27.1s 
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9.5 
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Abell 4059 
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0.015 


0.7 
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11980 
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-23d22ml0.4s 


0.0285 
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0.0 


28.8 


11980 


WBL 360-03 ^ 


Illi49m35.4s 
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0.0274 
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(1): Cluster name, (2-4): NED RA, Dec, redshift of BCG {htt-p: //nedwww. ipac. caltech. edu), (5): Reddening due to Galactic ex- 
tinction from Schlegel et al. (1998), (6): Spectroscopically- determined X-ray cooling rates (Mq yr~^ ) from McDonald et al. (2010), (7): 
1.4 GHz radio flux (mJy) from NVSS (http: // www. cv. nrao. edu/ nvss/ ) (8) HST proposal number for FUV data. Proposal Pis are W. 
Jaffe (#11131), W. Sparks (#11681), S. VeiUeux (#11980). 
^: No available Chandra data. 

O'Dea et al. 2008; hereafter MIR), and sub- mm (e.g., 
Edge 2001; Salome and Combes 2003) have allowed a 
much more detailed picture of star formation in BCGs. 
The typical star formation rates of ~ 1-10 Mq yr~^ 
(O'Dea et al. 2008) imply that gas at temperatures of 
~ 10^~^ K is being continuously converted into stars with 
an efficiency on the order of ~ 10%. The fact that most 
of these studies consider the integrated SF rates makes 
it difficult to determine the exact role of young stars in 
ionizing the extended warm gas observed at Ha, since 
the two may not be spatially coincident or the measure- 
ments may be contaminated by the inclusion of a central 
AGN. 

In order to understand both the role of star formation 
in ionizing the warm gas and the efficiency with which 
the cooling ICM is converted into stars, we have con- 
ducted a high spatial resolution far-UV survey of BCGs 
in cooling and non-cooling clusters. We describe the col- 
lection and analysis of the data from this survey in §2. 
In §3 we desribe the results of this survey, while in §4 
we discuss the implications of these results in the con- 
text of our previous work (M+10,M+11). Finally, in §5 
we summarize our findings and discuss any outstanding 
questions. Throughout this paper, we assume the follow- 
ing cosmological values: Hq = 73 km s~^ Mpc~^, O^rnatter 

= 0.27, ^vacuum — 0.73. 

2. DATA COLLECTION AND ANALYSIS 

To study the far-UV (hereafter FUV) emission in clus- 
ter cores, we selected 15 galaxy clusters from the larger 
samples of M+10 and M+11, which have deep, high spa- 
tial resolution (FWHM ~ 0.6^') Ha imaging from the 
Maryland-Magellan Tunable Filter (hereafter MMTF; 
Veilleux et al. 2010) on the Baade 6.5-m telescope at 
Las Campanas Observatory. Additionally, most of these 
systems have deep Chandra X-ray Observatory (here- 
after CXO) spectroscopic imaging, as well as Two-Micron 
Ah Sky Survey (hereafter 2MASS; Skrutskie et al. 2006) 
and NRAO VLA Sky Survey (hereafter NVSS; Condon 



et al. 1998) fluxes. This broad energy coverage provides 
an excellent complement to a FUV survey, allowing for 
the source of emission to be carefully identified. A sum- 
mary of these 15 clusters can be found in Table 1. For 
further information about the reduction and analysis of 
the Ha and X-ray data, see M+10. 

FUV imaging was acquired using the Advanced 
Camera for Surveys Solar Blind Channel (hereafter 
ACS/SBC) on the Hubble Space Telescope (hereafter 
HST) in both the F140LP and F150LP bandpasses when- 
ever possible, with a total exposure time of ^ 1200s 
each (PID #11980, PI Veilleux). The pointings were 
chosen, based on the results of our MMTF survey, to 
include all of the Ha emission in the field of view. Ex- 
posures with multiple filters are required to properly re- 
move the known ACS/SBC red leak, which has a non- 
negligible contribution due to the fact that the under- 
lying BCG is very luminous and red. Since the afore- 
mentioned filters are long-pass filters, they have nearly 
identical throughputs at longer wavelengths. Thus, by 
subtracting the F150LP exposure from the F140LP ex- 
posure we can effectively remove the red leak and con- 
sider only a relatively small range in wavelength, from 
1400A-1500A. Due to the small bandpass, it is possible 
for line emission to dominate the observed flux - we in- 
vestigate this possibility in §4. We have carried out this 
subtraction for 13/15 of the BCGs in our sample which 
have both F140LP and F150LP imaging. For Abeh 1795 
and Abell 2597, we are unable to remove the red leak 
contribution due to the lack of paired exposures, but we 
point out that, conveniently, these two systems have the 
brightest FUV flux in our sample and, thus, are largely 
unaffected by the inclusion of a small amount of non- 
FUV flux. Following the red-leak subtraction, we also 
bin the images 8x8 and smooth the images with a 1.5 
pixel smoothing radius, yielding matching spatial resolu- 
tion at FUV and Ha. This process is also necessary in 
order to increase the signal-to-noise of the FUV image, 
allowing us to identify interesting morphological features. 
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The final pixel scale for both the Ha and FUV images is 
0.2'7pixel. 

All FUV and Ha fluxes were corrected for Galactic 
extinction following Cardelli et al. (1989) using reddening 
estimates from Schlegel et al. (1998). 

3. RESULTS 

In the Appendix, we show the stellar continuum, Ha 
and FUV images for each of the 15 BCGs in our sam- 
ple. At a glance there does not appear to be consis- 
tent agreement between the Ha and FUV morpholo- 
gies. We observe systems having Ha filaments without 
accompanying FUV emission (Abell 1991, Abell 2052), 
systems with complex FUV emission without accompa- 
nying Ha (Abeh 1837, Abeh 2029), and systems with 
coincident Ha and FUV extended emission (Abell 1644, 
Abell 1795). Thus, it is obvious that a single explanation 
(e.g., star formation) is unable to account for the variety 
of FUV and Ha emission that we observe. 

As we did with the Ha emission in M+10, the FUV 
morphology can be classified as either nuclear or ex- 
tended. We find, in the FUV, 7/15 systems have ex- 
tended emission, 5/15 have nuclear emission, while 3/15 
have no emission at all. In order to quantitatively exam- 
ine both the nuclear and extended emission, we extract 
FUV and Ha fluxes in several regions, as shown in the 
Appendix. 

In Figure [l^, we show the correlation between the FUV 
and Ha luminosity for the regions identified in the Ap- 
pendix. We find a significant amount of FUV emission 
in all 5 of the systems for which we do not detect any 
Ha emission. Additionally, we see that at least 3 sys- 
tems are consistent with being shock-heated (Dopita & 
Sutherland 1996) - a point we will return to later in this 
section. 

As discussed by Hicks et al. (2010), a significant frac- 
tion of the FUV emission may be due to old stellar pop- 
ulations (e.g., horizontal branch stars) in the BCG. To 
proceed, we must isolate the FUV excess due to young, 
star- forming regions. In order to remove the contribu- 
tion from old stars, we consider the inner 3'^ and plot 
the K-band (from 2MASS; Skrutskie et al. 2006) ver- 
sus the FUV luminosity (Figure fT]3). Hicks et al. (2010) 
show that the FUV luminosity from old stars is highly 
concentrated in the central region, thus removing this 
contribution in the inner region will act as a suitable 
first-order correction. In order to calibrate this correc- 
tion for our sample, which lacks a control sample of con- 
firmed non-star- forming galaxies, we opt to fit a line 
which is chosen to pass through the four points with 
the lowest Lfuv/Lk' ratio. We make the assumption 
that these four galaxies with the lowest Lfuv/Lk' ra- 
tio are non-star-forming, which is supported by non- 
detections at Ha. The equation for this relation is: 
logio(LFt/y,3-) = 2.351ogio(LK',3-) - 42.98. The fact 
that four points with non-detections at Ha lie neatly 
along the same line suggests that this correction is mean- 
ingful. 

Figure [T]^ shows the FUV excess due to young stars 
versus the Ha luminosity for the total, nuclear and ex- 
tended regions in our complete sample of BCGs. With 
the contribution from old stellar populations removed, 
we find a tight correlation between Lpuv and Lna over 
four orders of magnitude. The majority of systems in 



our sample are consistent with the continuous star for- 
mation scenario (Kennicutt 1998, Leitherer et al. 1999), 
suggesting that much of the warm gas found in cluster 
cores may be photoionized by young stars. Two sys- 
tems, Abell 0970 and Abell 2029 have anomolously high 
FUV/Ha ratios, suggesting that star formation may be 
proceeding in bursts. As a starburst ages, the UV/Ha ra- 
tio will climb quickly due to the massive stars dying first. 
This means that, by 10 Myr after the burst, the UV/Ha 
ratio can already be an order of magnitude higher than 
the expected value for continuous star formation (see 
M+10 for further discussion). We find that the filaments 
in Abell 1991 and Abell 2052 are consistent with being 
heated by fast shocks, along with the nuclei of Abell 2052 
and Abell 2580. In the case of Abeh 2052, there ex- 
ist high quality radio and X-ray maps which show that 
the observed Ha emission is coincident with the inner 
edge of a radio-blow bubble. In Abell 1991, the Ha mor- 
phology is reminiscent of a bow shock, and is spatially 
coincident with a soft X-ray blob which is offset from the 
cluster core. Much of the FUV and Ha data is clustered 
between the regions depicting continuous star formation 
and shock heating, as shown in the zoomed-in portion of 
Figure [l] These regions may indeed be heated by a com- 
bination of processes, or they may simply be reddened 
due to intrinsic extinction. Based on their FUV/Ha ra- 
tios. Ha morphology, disrupted X-ray morphology, and 
high radio luminosity, we propose that the optical emis- 
sion in Abeh 1991, Abeh 2052, and Abeh 2580 is the 
product of shock heating, while the remaining 12 systems 
are experiencing continuous or burst-like star formation. 
We will return to this classification in the discussion. 
The SF rates that we measure in the systems with nu- 
clear emission only range from 0.01-0.1 yr~^, which 
are typical of normal red-sequence ellipticals (Kennicutt 
1998). For the systems with extended emission, exclud- 
ing those that are obviously shock-heated, the measured 
star formation rates range from 0.1-5 yr~^ which is 
similar to the rates of 0.008-3.6 Mq yr~^ observed in 
"blue early-type galaxies" (Wei et al. 2010). 

In Figure [2] we present the distribution of the FUV/Ha 
ratio in various regions for the 10 systems with detections 
(>l-<j) at either Ha or FUV. In the innermost region 
(r < 0.8^'), the warm gas appears to be shock- heated in 
60% of systems - these shocked nuclei may be associ- 
ated with AGN-driven outflows. Due to the small radial 
extent of this bin, the 2MASS data is of insufficient spa- 
tial resolution to remove any contribution from old stars. 
Thus, these FUV/Ha ratios are upper limits. Of the re- 
maining 4 systems, 2 are consistent with continuous star 
formation or a young starburst, while the remaining two 
are consistent with an aged starburst. At larger radius 
(0.8'' < r < 3.0'') the FUV/Ha ratio is slightly larger, 
with the distribution peaking in between the regions de- 
scribing shocks and star formation (see inset of Figure 

t. These data have had the contribution from old stel- 
r populations removed, as described earlier in this sec- 
tion. The FUV/Ha ratio at this radius is similar to what 
we observe in the filaments, as is seen in the third panel 
of Figure [2] The fact that the distribution of FUV/Ha 
peaks between the values for shocks and star formation 
supports a number of scenarios, including a mixture of 
the two processes, dusty star formation and star forma- 
tion with an IMF skewed towards high-mass stars (see 
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Fig. 1. — (a) FUV versus Ho; luminosity for 15 BCGs in our sample. For systems with extended FUV or Ho; emission, we separate the 
emission into extended, nuclear (inner 3^^) and entire systems (see Appendix). The dashed line represents the relation between Ho; and FUV 
luminosity for SF regions, as defined by Kennicutt (1998), the shaded blue region represents the expected FUV/Ho; ratio for continuous SF 
covering the full range of IMFs and metallicities from Starburst99 (Leitherer et al. 1999), while the shaded red region defines the expected 
FUV/Hq; ratio for fast shocks (Dopita &; Sutherland 1996). The legend in the lower right describes the different point type/colors, (b) 
FUV versus K^-band luminosity in the central ?>" of the BCG. The solid line represents our estimate of the contribution to the FUV 
luminosity from old, horizontal-branch stars, (c) Similar to panel (a), but with the contribution from old stars removed from the nuclear 
and total regions. The correlation between FUV and Ho; luminosity is much more significant, suggesting that the majority of the observed 
Hq; emission may be due to photoionization by young stars. The extended and nuclear emission in Abell 2052 and Abell 2580, as well as 
the extended emission in Abell 1991 appears to be consistent with heating by shocks, while the remaining systems are likely star-forming. 
The arrow in the lower right corner represents the magnitude of the intrinsic extinction correction for E(B-V)=0.2. 



right-most panel of Figure |2| . We will investigate these 
scenarios in §4. 

In M+10 and M+11, we showed that the Ha emission 
observed in the cool cores of galaxy clusters is intimately 
linked to the cooling ICM. In general, the thin, extended 
filaments observed in many of these clusters are found in 
regions where the ICM is cooling most rapidly, suggesting 
that this warm gas may be a byproduct of the ongoing 
cooling. If this is the case, it is relevant to ask what 
fraction of the cooling ICM is turning into stars. We 
address this question in Figure [S] by comparing the star 



formation rate with the X-ray cooling rate (dM/dt) for 32 
galaxy clusters. In order to compute the star formation 
rate, we use the prescriptions in Kennicutt (1998). For 
the systems observed with HST, we use the average of 
the FUV- and Ha-determined star formation rates (filled 
blue circles). For an additional 10 clusters from M+10 
and M+11 we make use of archival GALEX data for 5 
clusters (open blue circles) and assume that both the UV 
and Ha emisison trace star forming regions. These data 
have also had the contribution from old stellar popula- 
tions removed, as described in M+10. In the remaining 
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Fig. 2. — Distribution of the FUV to Ho; luminosity ratios for 
the 10 clusters in our sample with high S/N detections in either 
the FUV or Ho;. Each panel considers a different region, from 
top to bottom: the nucleus, the central 3" , excluding the nucleus, 
the filaments, and the entire system. The red vertical band rep- 
resents the region consistent with shocks (Dopita &: Sutherland 
1996), while the blue band represents the region consistent with 
continuous star formation (Leitherer et al. 1999). In the nucleus, 
the majority of the observed Ho; emission is consistent with having 
been shock-heated, while in the outer regions, including filaments, 
the FUV/Hq; ratio lies between the two highlighted regions. In the 
bottom panel, the dotted histogram shows how the distribution 
would be altered if we corrected the fluxes for an i intrinsic ex- 
tinction of E(B-V)=0.2 (dust screen model), while the blue dashed 
line represents the expected value for 08V stars (symbolic of a 
top-heavy IMF). 



5 cases where there is no accompanying UV data, we as- 
sume that the Ha emission is the result of photoioniza- 
tion by young stars, and convert the Ra luminosity into 
a continuous star formation rate (green triangles). For 
the three shock-heated systems (red crosses) mentioned 
above, we determine the SF rate based on the FUV lu- 
minosity alone - this represents an upper limit on the 
amount of continuous star formation. Finally, we also in- 
clude 10 clusters from the sample of O'Dea et al. (2008), 
who compute SF rates based on Spitzer dditdi, and dM/dt 
in a similar manner to us (open purple circles). We note 
that the SFRs derived from FUV data are systematically 
lower than those derived from 24/im Spitzer data. This 
is most likely due to our lack of an intrinsic extinction 
correction for these FUV data. 

We find that the "efficiency" of star formation, defined 
as the current ratio of stars formed to gas cooling out 
of the ICM, can range from 1% to 50% over the full 
sample. This spread is independent of whether we only 
consider systems classified as star-forming based on their 
FUV/Ha ratios in Figure [l] For the majority of the 32 
clusters in Figure js) there is a tight correlation between 
SF rate and dM/dt with a typical efficiency between 10- 
50%. In the right panel of Figure [sj we show that the 
efficiency is nearly independent of the cluster tempera- 
ture, with only a weak dependence which is primarily 
driven by highly extincted systems. This suggests that 
the temperature of the surrounding ICM does not ham- 
per the BCG's ability to form stars. While not shown 
here, we also investigated the distribution of star for- 
mation efficiencies with the central entropy, Kq^ from 
Cavagnolo et al. (2009) and, similarly to Tx, ffnd no cor- 
relation. In the following section we will discuss possible 
interpretations of this efficiency measure. 

4. DISCUSSION 

4.1. Star Formation as an Ionization Source 

In Crawford et al. (2005), a variety of ionization 
sources are discussed which could produce the observed 
Ha emission in the cool cores of galaxy clusters. The 
purpose of this HST survey was to investigate one of 
the most plausible ionization sources: photoionization 
by young stars. In Figures [l] and [2] we showed that, 
once the contribution to the FUV emission from old stel- 
lar populations is removed, the majority of the Ha and 
FUV emission that we observe in cluster cores is roughly 
consistent with the star formation scenario. Based on 
the FUV/Ha ratios, we identify three different types of 
system: 

1. FUV/Ha > 10-^2 Hz-^ Suggests a starburst 
that has aged by at least 10 Myr. Two systems, 
Abeh 0970 and Abeh 2029, fulffh this criteria, while 
several others may fall into this category if their Ha 
luminosity is significantly less than the measured 
upper limits. 

2. FUV/Ha - 10-13 Rz'^. The FUV/Ha ratios of 
these systems are consistent with continuous star 
formation or a recent (0-5 Gyr ago) burst of star 
formation. The filaments in Abell 1644, Abell 1795, 
Abeh 2597, and part of Abeh 2052, along with the 
nuclei of Abeh 1795, Abeh 1991, Abeh 2151, and 
Abell 2597 appear to be star-forming. 
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Fig. 3. — Left: SF rate (SFR) versus X-ray cooling rate (dM/dt) for 32 galaxy clusters. The filled blue circles (10) and red crosses (2) 
correspond to clusters identified as star-forming and shock-heated from Figure^ respectively. One of the shock-heated systems, Abell 2580, 
does not have available data from the CXO archives. The open blue circles (5) refer to clusters with MMTF Ho; imaging and archival 
GALEX near-UV (NUV) data from M+10. The green triangles (5) are clusters with deep Ho; imaging from M+10, M+11 - we assume this 
emission is due to star formation. Star formation rates are computed from NUV, FUV, and Ho; luminosities based on the prescriptions in 
Kennicutt (1998). We also show data from O'Dea et al. (2008) as open purple circles (10), with SF rates based on mid-IR (MIR, 24/im) 
Spitzer data and dM/dt computed in a similar fashion to us. The blue star represents the Perseus cluster (Conselice et al. 2001, Sanders 
et al. 2004). Diagonal lines refer to efficiencies of 1%, 10% and 100%, while the shaded regions refer to the typical SF efficiency in molecular 
clouds per free fall time (gray; Krumholz Sz McKee 2005) and over the lifetime of the cloud (tan; Kroupa et al. 2001; Lada & Lada 2003). 
Typical errors (20% in SFR, 50% in dM/dt) are shown in the bottom right corner. Right: SF efficiency (SFR per dM/dt) as a function of 
cluster X-ray temperature from the ACCEPT sample of Cavagnolo et al. (2009). Typical errors are shown in the upper right corner. We 
note that the error in SF efficiency is dominated by our uncertainty in the X-ray cooling rate, dM/dt. The highest-temperature point is 
the Ophiuchus cluster, which lies close to the Galactic plane and, thus, suffers from heavy extinction. While we do correct for Galactic 
extinction, the large column density will amplify any uncertainty in the measurement of Nh or in the extinction curves applied. 



3. FUV/Ha < 10"^"^ Rz'^. Suggests heating by fast 
shocks or some other source of hard ionization (e.g., 
cosmic rays, AGN). The filaments of Abeh 2052 
and Abeh 1991, and the nuclei of Abeh 2052 and 
Abell 2580 have FUV/Ha ratios which are consis- 
tent with this picture, in the absence of internal 
reddening. 

Figures [l] (inset) and [2] show that a large fraction of 
the systems which we oDserve fall between the regions 
describing shock heating and star formation. However, 
these data have not been corrected for intrinsic redden- 
ing due to dust. Correcting for a very modest reddening 
(E(B-V) - 0.2) would boost the FUV luminosity of these 
systems such that the FUV/Ha ratio is consistent with 
star formation (see Figure ITp and the lower panel of Fig- 
ure |2|. Unfortunately, the amount of reddening in the 
filaments and nuclei of these systems is currently not well 
constrained for very many systems, but typical values of 
E(B-V) can range from 0-0.4 in the cores of galaxy clus- 
ters (Crawford et al. 1999). In the case of Abeh 2052, 
for which we measure FUV/Ha ratios indicative of shock- 
heating and have an estimate of the amount of intrinsic 
reddening from Crawford {E{B — V) = 0.221q'2^), we can 
investigate whether correcting for this extinction would 
provide FUV/Ha ratios consistent with star- forming re- 
gions. Assuming a simple dust-screen model, correcting 
for a reddening of E{B — V) = 0.22 would transform a 



FUV/Ha ratio of 4.4 x 10"^^i^2:"^ in the filaments of 
Abeh 2052 to 1.2 x 10"^'^i^^"\ which is consistent with 
the upper limit for shock-heated systems (see Figure [2|. 
However if, contrary to expectations, the filaments have a 
slightly higher reddening than the nucleus, the FUV/Ha 
ratio may be even higher. Thus, it is certainly possi- 
ble that the systems which we classify as shock-heated, 
or those which have ambiguous FUV/Ha ratios, may in 
fact be highly-obscurred star- forming systems. We will 
address this possibility in significantly more detail in an 
upcoming paper which will include long-slit spectroscopy 
of the Ha filaments providing, for the first time, redden- 
ing estimates away from the nucleus in these systems. 

An alternative explanation for the intermediate 
FUV/Ha ratios is that the IMF in the filaments is top- 
heavy {a <C 2.35). Again, the lower panel of Figure [2] 
shows that the peak of the FUV/Ha distribution is con- 
sistent with the value expected for 08V stars. There is 
a substantial amount of literature providing evidence for 
a top-heavy IMF in various environments including the 
Galactic center (Maness et al. 2007) and disturbed galax- 
ies (Habergham et al. 2010). Thus, regardless of whether 
there is a small amount of dust or a slightly altered IMF, 
we suspect that the majority of the systems with inter- 
mediate FUV/Ha ratios are in fact star- forming, with 
the exception of Abeh 1991, Abeh 2052, and Abeh 2580, 
which have low FUV/Ha ratios and morphologies which 
resemble bow-shocks and/or jets. 
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Fig. 4. — Left panel: UV spectral energy distributions (SEDs) for 6 cluster cores in our sample using data from GALEX and XMM-OM 
(M+10). All fluxes have been corrected for Galactic extinction. The two GALEX fllters are denoted by G"*^ and G"^, while the three XMM- 
OM fllters are denoted by X-*^, X^ and X^. The colored lines represent a flt to the GALEX+XMM-OM data, while the stars represent the 
FUV flux that we measure with HST ACS/SBG. Right panel: Residuals from the continuum flts shown in the left panel. Our new HST 
FUV data are consistent with the UV continuum flt, showing little evidence for additional line emission. 



Due to our use of the F150LP filter to remove red 
leak contamination, we are considering only a very small 
wavelength range from 1400-1500A. In this region, there 
may be line emission from [OIV] and various ionization 
states of sulfur due to gas cooling at ~ 10^ K. In or- 
der to establish that we are indeed observing continuum 
emission from young stars, we have computed UV spec- 
tral energy distributions (SEDs) for 6 BCGs which have 
deep GALEX, XMM-OM, and HST UV data. These 
data are presented in the left panel of Figure [4j We see 
that, in general, the UV SED follows a powerlaw over the 
range of 1500-3000A. The new HST data, depicted as 
colored stars in this plot, agree well with the extrapola- 
tion of the continuum to shorter wavelengths, suggesting 
that there is very little contamination from line emission. 
This also suggests that there is little contribution from 
a diffuse UV component. This is further emphasized in 
the right panel of Figure [4] where we show the residu- 
als from the continuum fit lor each of the 6 BCGs. Our 
measured FUV fluxes from these new HST data are con- 
sistent with the measurement of a UV continuum from 
archival GALEX and XMM-OM data. 

The idea that massive, young stars may be responsible 
for heating the majority of the warm, ionized filaments 
observed in cool core clusters is certainly not a new one 
(see e.g., Hu et al. 1985; Heckman et al. 1989; McNamara 
and O'Connel 1989). Most recently, O'Dea et al. (2008), 
Hicks et al. (2010) and McDonald et al. (2010a) con- 
ducted MIR, UV and Ha surveys, respectively, of cool 
core clusters and found a strong correlation between 
the SF rate and the cooling properties of cluster cores. 
However, this work extends these findings to include 
spatially-resolved SF rates, which the previous studies 
have been unable to provide. This allows us to say con- 
clusively that the young stars and the warm, ionized gas 



are in close proximity (< f ) in the vast majority of 
systems, offering a straightforward explanation for the 
heating of these filaments. 

4.2. Star Formation Efficiencies in Cooling Flows 

In §3, we provide estimates of the efficiency with which 
the cooling ICM is converted into stars, assuming that 
this is indeed the source of star formation. This as- 
sumption is based on the results of M+10 and M+11, 
which provided several strong links between the X-ray 
cooling properties and the warm, ionized gas. These es- 
timates of star formation efficiency represent a constraint 
on the so-called "accreting box model" of star formation. 
The simplified model that we propose is that the ICM 
is allowed to cool rapidly in regions where cooling lo- 
cally dominates over feedback (Sharma et al. 2010). Our 
estimates of the ICM cooling rate, based on medium- 
resolution CXO spectra, are consistent with estimates 
based on high-resolution XMM grating spectroscopy by 
Peterson et al. (2003) for the five overlapping systems. 
Once the gas reaches temperatures of ~10^"^K, it can 
continue to cool rapidly via UV/optical/IR line emis- 
sion without producing fluxes that are inconsistent with 
what are observed. In the standard way, star formation 
will proceed once the gas reaches low enough temper- 
ature and high enough density. Observations by Edge 
(2001) and Salome & Combes (2003) show evidence for 
molecular gas in the cool cores of several galaxy clusters, 
consistent with this picture. 

If the above scenario is correct, our estimate of the 
star formation efficiency provides a constraint on the 
fraction of hot gas that will be converted to stars, as- 
suming a steady inflow of gas. In Figure [5] we provide a 
histogram of SF efficiencies for all of the systems in Fig- 
ure 3] with non-zero X-ray cooling rates (dM/dt). The 
peakof this distribution is well defined at an efficiency of 
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Utl %, regardless of which SF indicator (UV, Ha, MIR) 
is used. We note that, while the distribution for UV- and 
MIR-determined SF rates both peak at roughly the same 
value, the UV-determined SF histogram extends to much 
lower values. The low-efficiency tail of this distribution 
may be an artifact produced by intrinsic extinction due 
to dust, to which the UV will be most sensitive, or may 
be indicative of a selection bias in the MIR sample. If 
we measure the peak efficiency based on the subsamples 
excluding the MIR and MIR+Ha (with no accompany- 
ing UV) data, we get 10tf% and 14tf %, respectively. 
Thus, the peak value of 14% is not solely driven by the 
inclusion of MIR data. 

The average efficiency of 14^^^%, based on MIR, Ra 
and UV data, is consistent with the estimates of star 
formation efficiency over the lifetime of a typical molec- 
ular cloud (20-50%; Kroupa et al. 2001, Lada & Lada 
2003). This large variance in star formation efficiency 
may be due to differences in the ICM cooling and star 
formation timescales. Naturally, one would expect that 
there is some delay between the ICM cooling and the 
formation of stars, so that a reservoir of cold gas can 
accumulate and the formation of stars can be triggered. 
If this is indeed the case, one would expect to observe 
cooling-dominated periods (low SFE) followed by peri- 
ods of strong star formation (high SFE) once the cold gas 
reservoir has reached some critical mass. Over an ensem- 
ble of systems, the average SFE is then an estimate of the 
time-averaged efficiency of an accreting system in con- 
verting a steady stream of cooling gas into stars. An al- 
ternative explanation for the spread of observed efficien- 
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Fig. 5. — Distribution of measured star formation efficiencies 
from Figure Is] for 26 systems. Systems which are hkely shock- 
heated (Abeiri991, Abeh 2052, Abeh 2580) have not been in- 
cluded. The additive contribution to the total histogram (solid 
black line) from UV+Ho;, Ho; and MIR data are shown in differ- 
ent colors following the color scheme in Figure [3] The dotted line 
shows a Gaussian fit to this histogram, which peaks at 14^g^% 
efficiency. 



cies is that the source of feedback is episodic (e.g., AGN). 
In this scenario, an episode of strong feedback from the 
AGN would re-heat the reservoir of cool gas, severely 
reducing the potential for star formation. This may in- 
deed be the case, since two of the three systems with the 
highest 1.4 GHz luminosity (Abeh 2052, Abeh 2597, and 
Perseus A) have SFE < 0.1. 

The fact that Figure [5] shows a well-defined peak sug- 
gests that the fraction of stars formed in an accreting 
system is constant over long enough timescales. Our 
estimate of an average efficiency indicates that, for a 
steady-state system accreting hot gas which is then al- 
lowed to cool, roughly 4 Mq of gas will either be re- 
heated or expelled via winds for every 1 Mq of stars 
formed. This fraction of baryons in stars is consistent 
with the global fraction of ~ 20-30% required by simula- 
tions to reproduce the observed stellar mass function of 
galaxies (Somerville et al. 2008). Unlike measurements 
of SF efficiency for giant molecular clouds, this estimate 
does not require the use of a specific timescale, since we 
are assuming that stars are forming out of the inffow of 
hot gas and that the reservoir for this hot gas is inex- 
haustible. 

5. SUMMARY AND FUTURE PROSPECTS 

We have assembled a unique set of high spatial res- 
olution far-UV and Ra images for 15 cool core galaxy 
clusters. These data provide an unprecedented view of 
the thin, extended fflaments in the cores of galaxy clus- 
ters. Based on the ratio of the far-UV to Ha luminosity, 
the UV SED, and the far-UV and Ha morphology, we 
conclude that the warm, ionized gas in the cluster cores 
is photoionized by massive, young stars in all but a few 
(Abeh 1991, Abeh 2052, Abeh 2580) systems. We show 
that the extended filaments, when considered separately, 
appear to be forming stars in the majority of cases, 
while the nuclei tend to have slightly lower FUV/Ha 
ratios, suggesting either a harder ionization source or 
higher extinction. The slight deviation from expected 
FUV/Ha ratios for continuous star formation (Leitherer 
et al. 1999) may be due to the fact that we have made no 
attempt to correct for intrinsic extinction due to dust or 
due to a top-heavy {a <C 2.35) IMF. We note that mod- 
est amounts of dust (E(B-V) ~ 0.2) in the most dense 
regions of the ICM can account for this deviation. Ide- 
ally, one would like spatially-resolved optical spectra of 
the filaments in order to constrain the heat source and 
intrinsic reddening of the filaments. We plan on address- 
ing this issue in upcoming studies. Comparing the esti- 
mates of the star formation rates based on FUV, Ha and 
MIR luminosities to the spectroscopically-determined X- 
ray cooling rate suggests a star formation efficiency of 
14^3^%. This value represents the time- averaged frac- 
tion, by mass, of gas cooling out of the ICM which turns 
into stars and agrees well with the stars-to-gas fraction 
of ~20-30% required by simulations to reproduce the ob- 
served stellar mass function. This result provides a new 
constraint for studies of star formation in accreting sys- 
tems. Many aspects of this simplified scenario are still 
not well understood, including whether the star forma- 
tion is similar to that seen in nearby spirals or vastly 
different. We intend to investigate such differences via 
an assortment of star formation indicators from the UV 
to radio in future work. 
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APPENDIX 




Fig. 6. — Optical and FUV data for the 15 clusters in our sample. From left to right the panels are: 1) MMTF red continuum image, 
2) MMTF continuum-subtracted Ra image, 3) ACS/SBC FUV redleak-subtracted image (F150LP-F140LP), 4) ¥La image with extraction 
regions defined. The horizontal scale bar in the left two panels represents 20 kpc. The Ha and FUV images are zoomed in relative to the 
red continuum image. The square region in the red continuum panels represents the field of view for the zoomed-in panels. The grayscale 
in all images is arbitrarily chosen in order to enhance any morphological features. 




Fig. 



6. — Continued. 




Fig. 6. — Continued. 



